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ABSTRACT 


In  this  paper,  the  authors  reviewed  the  literature  on  computa¬ 
tional  aspects  of  the  distributions  of  the  likelihood  ratio  statistics 
for  testing  various  hypotheses  on  the  covariance  matrices  and  mean 
vectors  of  complex  multivariate  normal  populations.  Applications  of 
some  of  these  test  procedures  in  the  area  of  inference  on  multiple 
time  series  in  the  frequency  domain  are  also  discussed.  In  the  Appendix, 
the  authors  give  tables  which  are  useful  in  implementation  of  various 
likelihood  ratio  test  statistics  discussed  in  this  paper. 

Key  words  and  phrases; 

Complex  multivariate  normal,  homogeneity  of  populations,  likelihood 
ratio  tests,  multiple  independence,  multiple  time  series  in  the  frequency 
domain,  and  tables. 
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1 .  Introduction 


It  is  known  that  certain  estimates  of  spectral  density  matrices  of 
the  stationary  and  Gaussian  multiple  time  series  are  distributed  approxi¬ 
mately  as  complex  Wlshart  matrices.  So,  complex  multivariate  distributions 
are  useful  (e.g.,  see  Brlllinger  (1974),  and  Hannan  (1970))  in  the  area 
of  inference  on  multiple  time  series.  These  distributions  are  useful  in 
the  area  of  nuclear  physics  (see  Carmell  (1974))  also. 

Wooding  (1956)  Introduced  the  complex  multivariate  normal  distribution. 

A  complex  random  vector  Is  said  to  be  distributed  as  a  complex  multivariate 
normal  if  its  real  and  imaginary  parts  are  distributed  jointly  as  a  multi¬ 
variate  normal  with  a  structered  covariance  matrix.  Motivated  by  applica¬ 
tions  In  time  series,  Goodman  (1963a,  b)  made  a  systematic  study  of  the 
complex  multivariate  normal  distribution  and  complex  Wlshart  matrix.  Since 
then,  James  (1964),  Khatrl  (1965),  Krlshnaiah  (1976)  and  other  workers  in 
the  field  have  investigated  various  aspects  of  complex  multivariate  distribu¬ 
tions.  For  a  review  of  the  literature  on  complex  multivariate  distributions, 
the  reader  is  refered  to  Krlshnaiah  (1976).  In  this  paper,  we  review  the 
literature  on  the  likelihood  ratio  tests  on  mean  vectors  and  covariance 
matrices  of  the  complex  multivariate  normal  populations  as  well  as  some  of 
their  applications  in  the  area  of  Inference  on  multiple  time  series  in 
the  frequency  domain. 

In  Section  2  of  this  paper,  we  discuss  the  complex  multivariate  normal 
and  complex  Wlshart  matrix.  The  distribution  of  the  determinant  of  the  com¬ 
plex  multivariate  beta  matrix  is  discussed  in  Section  3,  whereas  Section  4 
is  devoted  to  the  likelihood  ratio  test  procedure  for  testing  the  hypothesis 
of  multiple  Independence  of  several  sets  of  variables  when  their  joint  distribu¬ 
tion  is  complex  multivariate  normal.  Likelihood  ratio  tests  for  the  hypothesis  of 
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sphericity  and  the  hypothesis  specifying  the  covariance  matrix  are  discussed 
In  Sections  5  and  6  respectively.  In  Section  7  we  discuss  the  likelihood 
ratio  test  for  the  homogeneity  of  the  covariance  matrices  whereas  the 
likelihood  ratio  test  procedure  for  the  homogeneity  of  several  complex 
multivariate  normal  populations  Is  discussed  In  Section  8.  Likelihood 
ratio  test  procedure  specifying  the  covariance  matrix  and  mean  vector  Is 
discussed  In  Section  9.  Applications  of  some  test  procedures  on  the  covari¬ 
ance  matrices  of  the  complex  multivariate  normal  populations  to  the  area  of 
Inference  on  multiple  time  series  in  the  frequency  domain  are  discussed  In 
Section  10.  Various  tables  useful  In  Implementation  of  certain  likelihood 
ratio  test  procedures  are  given  in  the  Appendix.  These  tables  are  constructed 
by  approximating  a  suitable  power  of  the  likelihood  ratio  statistics  with 
Pearson's  Type  I  distribution  by  using  the  first  four  moments.  The  accuracy 
of  these  tables  Is  found  to  be  good. 
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2.  Complex  Multivariate  Normal  and  Complex  Wlshart  Distributions 


Let  z  •  X  +  ly  where  x  and  y  are  of  order  p^l  and  (x',  y')  is  distri¬ 
buted  as  2p-variate  normal  with  mean  vector  and  covariance  matrix 

c.h 

where  A'  denotes  the  transpose  of  A.  Then,  the  distribution  of  z  is  known 
to  be  the  complex  multivariate  normal  distribution  with  mean  vector  p  and 
covariance  matrix  E  where  u  +  E“2(E^-  112)*  probability 

density  function  (p.d.f.)  of  z  is  given  by 

f(z)  =  — ^ —  exp{  -(z-p)  'r  ^  (z-p)}  (2.2) 

ttPIzI 

whereas  the  characteristic  function  of  z  is 

«(t)  -=  expdt'p  -  t'r  (2.3) 

where  t  =  tj^+it2  and  t  is  the  complex  conjugate  of  t.  Wooding  (1956)  derived 
expressions  for  the  p.d.f.  and  characteristic  functior  of  z.  The 


(2.1) 


Maximum  likelihood  estimates  of  p  and  E  based  on  a  random  sample 
(zj,...z^)  are  known  to  be 


N 


N-1  I 


E  •  n“^  ^(z,-p)  (z.-p) ' 


N 


-j 


j-1 


tj  r'  'tj 


(2.4) 
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Also,  V  and  Z  are  distributed  independent  of  each  other. 


Next,  let  S  =  NZ.  Then,  the  distribution  of  S  is  known  to  be  a  central 
complex  Wlshart  matrix  with  n»N-l  degrees  of  freedom.  The  probability  density 
of  S  is  known  (Goodman  (1963b))  to  be 


f(S) 


in-p 


etr{-Z 


tP(p"1)/2  n  r(n-j+l)  lz|” 


(2.5) 


j-1 


where  ^tr  B  denotes  the  exponential  of  the  trace  of  the  matrix  B. 
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3.  Distribution  of  the  Determinant  of  the  Complex  Multivariate  Beta 
Matrix 


In  this  section,  we  discuss  the  distribution  of  the  determinant  of  the 
complex  multivariate  beta  matrix.  This  distribution  in  useful  for  testing 
the  hypothesis  of  the  equality  of  several  mean  vectors  and  the  equality  of 
two  covariance  matrices  when  the  underlying  distributions  are  complex  multi¬ 
variate  normal.  It  is  also  useful  in  testing  the  hypothesis 

-  —  - 

where  covariance  between  two  sets  of  variables  whose  joint  distri¬ 

bution  is  complex  multivariate  normal. 

Let  A^s  P P  and  A25  p x p  be  independently  distributed  as  the  central  com¬ 
plex  Wishart  matrices  with  n  and  q  degrees  of  freedom,  and  let  E(Aj^/n)  = 

E(A2/q)  =  I.  Then  Aj^(Aj^+A2)  ^  is  known  to  be  a  (central)  complex  multivariate 
beta  matrix.  Now,  let 

U  -  |a^(A2  +  A2)~^|.  (3.1) 


The  h-th  moment  of  U  is  known  to  be 


E(U^) 


P 

n 


r(n+h-1+l)  r(n-fq-1+l) 
r(n-j+l)  r(n+h4q-j+l) 


(3.2) 


Using  the  first  four  moments  of  Lee,  Krishnaiah  and  Chang  (1975)  have 

approximated  the  distribution  of  with  the  Pearson's  Type  I  distribution 

where  b  is  a  properly  chosen  Integer.  The  constant  b  is  chosen  to  be  equal 
to  1  or  2  according  as  M >  20  or  M <  20  where  M  »  n  -  p+1 .  By  making  use 
of  this  approximation,  values  of  Cj^  are  computed,  where 


PlCi^Cil  «  (1-a), 


(3.3) 


z  z 

C^  *  -  (2n  +  q-p)  log  U/X2pq  and  X2pq  ^  is  the  upper  100  a  /o  value  of 
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X  with  2pq  degrees  of  freedoo.  The  values  of  are  computed  for  a  -  0.005, 
0.01.  0.025,  0.05,  0.1,  0.90,  0.95.  0.99,  0.995,  M  -  1(1)10(2)20,  30,  60,  120, 
where  M  ■  n-p-fl.  These  values  are  given  In  a  technical  report  by  Lee, 
Krlshnalah  and  (^ng  (1975).  The  upper  5Z  and  IZ  points  are  reproduced  in 
Table  7  of  this  chapter.  To  check  for  the  accruacy  of  the  entries  in  Table 
7 ,  the  above  authors  compared  some  of  the  values  obtained  by  the  Pearson  type 
approximation  with  the  corresponding  exact  values.  These  comparisons  are  given 
in  Table  1. 


Table  1 

Comparison  of  the  Pearson  type  Approximation  with  exact  expression  for 
the  Distribution  of  c^ 


P 

■  2  q  * 

3 

P 

■2  q  ■  i 

20 

M 

a 

L-K-C 

Exact 

a 

L-K-C 

Exact 

1 

0.05 

1.286 

1.289 

0.05 

1.928 

1.932 

1 

0.01 

1.350 

1.349 

0.01 

2.085 

•2.080 

5 

0.05 

1.029 

1.029 

0.05 

1.243 

1.243 

5 

0.01 

1.033 

1.033 

0.01 

1.262 

1.262 

9 

0.05 

1.010 

1.011 

0.05 

1.12? 

1.128 

9 

0.01 

1.012 

1.012 

0.01 

1.137 

1.137 

The  constant  a  in  Table  1  is  defined  by  Eq.  (3.3).  Also,  the  values  under 
the  column  "L-K-C"  are  the  values  of  c^  obtained  by  Lee,  Krishnaiah  and 
Chang  (1975)  using  the  Pearson  type  approximation  whereas  the  values  under 
the  column  "Exact"  are  the  corresponding  values  given  by  Gupta  (1971).  Table 
1  indicates  that  the  accruacy  of  the  Pearson  type  approximation  is  sufficient 
for  practical  purposes. 
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4.  Test  for  Independence  of  Sets  of  Variates 

I  '  ,  ^ 

Let  z  »  (Zj  )  be  distributed  as  a  complex  multivariate  normal 

•  '  t 

distribution  with  mean  vector  P  »  )  and  covariance  matrix  1.  Also, 

r  _ 5  -  .1 

let  -  y^)  ^ij  ’  assumed  that  is  of 

order  p^ 1  and  p^+. . . +p^  »  s.  In  this  section  we  discuss  the  problem  of 
testing  the  hypothesis  H2  where 


“2=  ^Ij  ■  0 


(1.1) 


for  i?ij=l . q.  Now  let 


A  = 


All  Aj^2  ••• 

A21  A22  ...  A2q 


I  Am  •••  A 

L  qq 


where 


=  y  (i.-z  )(z,.-2,),2  =N  y 

^  '-ffi  .h  ’  o  ^ 


"gh  tg/'thj  ‘h.-  ’  -g. 


“gj 


and  (7 z'  j)  Is  j-th  independent  observation  on  (z’,...,z'). 
-Ij  ~  qj  -1  -q 

The  likelihood  ratio  statistic  for  testing  H2  is 

lAl 


q 

n  |A 


(4.2) 


ij 


In  the  analogous  real  case,  Wilks  (1935)  derived  the  likelihood  ratio  test 
for  multiple  independence.  The  hth  moment  of  >2 

qnjigjLl-).. 

j,!  r(n-j;i)  J  r(n+h-a+l) 


£(^2) 


(4.3) 
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where  n =  N-1  and  r(*)  is  the  complete  gamma  function.  The  distribution  of 
1/4 

^2  is  approximated  by  Pearson's  type  I  distribution  .with  density 

g(x)  »  {B(oi+l,e+l)  ^(x-aQ)“(aj^-x)^  (4.4) 

where  Oq  <  x  <  and  a  and  e  are  some  real  munbers. 

Approximate  percentage  points  of  the  distribution  of  =  -2  log  >2 
constructed  by  Krishnaiah,  Lee  and  Chang  (1975,  1976)  for  Pj^=p=l,2,3;  q=3,4,5» 
a  =  0.01,  0.05,  0.10;  M  =  1(1)20(2)30  where  M  *  n-s-3,and 

Pr[X2  £  C2IH2]  =  (1-a).  These  percentage  points  are  reproduced  in  Table  8. 
Percentage  points  for  q=2  are  given  in  Table  7. 

Now,  consider  a  class  of  statistics  W(0  £  W  <_  1)  whose  moments  are.  of 

..  .  a 

r(x^(l+h)  +4^] 

- ,  h=0,l,...  (4.5) 

n  r[y  (l+h)+y  ] 
j-1  ^  ^ 

a  c 

where  K  is  a  normalizing  constant  such  that  E{W°}  *  1  and  V  =  T  y.. 

k=l  j=l  ^ 

The  likelihood  ratio  test  statistics  considered  in  this  chapter  are 
special  cases  of  the  above  class  of  statistics.  Box  (1949)  gave  explicitly 
the  first  few  terms  of  an  asymptotic  expression  for  the  distribution  of  a 
class  of  statistics  whose  moments  are  of  the  form  (4.5).  But  the  first  few 
terms  alone  are  not  sufficient  to  get  the  desired  degree  of  accuracy  in 
a  number  of  practical  situations.  So,  Lee,  Krishnaiah  and  Chang  (1976) 
gave  terms  up  to  order  n  explicitly.  In  Table  2,  given  below,  a 
comparison  of  the  values  obtained  by  using  Pearson  type  approximation  is 

made  with  the  corresponding  values  obtained  by  the  asymptotic  expression 

c  j  -13 
of  order  n 


Efv/'}  =  K  n  y^^  n 
j-1  3 
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Table  2 


Comparison  of  the  Pearson  Type  Approximation  with  the  Asymptotic 
Expansion  for  the  Distribution  of 


n 

q  = 

3  P 

-  1 

q 

-  4 

P  -  2 

— 

q  “ 

5  P 

-  1 

^^2 

'  “1 

“2 

*^2 

“1 

“2 

*^2 

“1 

“2 

10 

1.459 

0.05 

0.0499 

4.011 

0.05 

0.0437 

10 

1.949 

0.01 

0.0100 

- 

- 

- 

4.811 

0.01 

0.0095 

15 

0.923 

0.05 

0.0500 

5.733 

0.05 

0.0479 

2.435 

0.05 

0.0497 

15 

1.233 

0.01 

0.0100 

6.496 

0.01 

0.0093 

2.914 

0.01 

0.0099 

20 

0.675 

0.05 

0.0500 

3.958 

0.05 

0.0493 

1.752 

0.05 

0.0499 

20 

0.902 

0.01 

0.0100 

4.479 

0.01 

0.0098 

2.096 

0.01 

0.0100 

30 

0.439 

0.05 

0.0501 

2.455 

0.05 

0.0498 

1.124 

0.05 

0.0500 

30 

0.587 

0.01 

0.0100 

2.777 

0.01 

0.0099 

1.344 

0.01 

0.0100 

In  this  table,  is  the  value  of  a  if  we  use  the  Pearson  type 

approximation  and  02  is  the  value  of  a  if  we  use  the  asymptotic  expression 
-13 

of  order  n  .  From  the  table,  we  observe  that  the  accuracy  of  the  Pearson 
type  approximation  is  sufficient  for  practical  purposes. 
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5.  Test  for  Sphericity 


The  likelihood  ratio  statistic  for  testing  :  I  =  o  is  given  by 

.-1| 


A  Z 


^3  = 


(tr  A  Z  ^/s)® 
o 


(5.1) 


where  Z  is  known,  A  was  defined  in  Section  ^  and  trA  denotes  the  trace  of 
o 

A.  The  h^^  moment  of  known  to  be 


„,.h,  r  (sn)  I  r(n  +  h-.i  +1)  .  (5.2) 

°  r(sn  +  sh)  r(n-j+l) 

Mauchly  (1940)  derived  the  likelihood  ratio  statistic  for  testing  the  hypothe 
sis  df  sphericity  when  the  underlying  distribution  is  real  multivariate  normal 

l/b 

The  distribution  of  is  approximated  by  a  Pearson  Type  I 
distribution,  where  b  is  a  suitably  chosen  integer.  For  M^22,  we  took 
b=2  and  for  M<22,  b=>4.  Using  the  approximation  described  above,  approxi¬ 
mate  upper  percentage  points  of  distribution  of  X^  =-2  log  X^  were 
constructed  by  Krishnaiah,  Lee  and  Chang  (1975,  1976)  for  s  =  2(1)10, 
a  =  0.01,  0.05,  M  =  1(1)  20(2)30(5)50,60,  where  M  =  n-s-3  and 
P[X3  _<  C3IH3]  =  (1  -  a).  These  values  are  reproduced  in  Table  9. 

In  Table  3,  we  compare  the  values  obtained  by  the  Pearson  type  approx¬ 
imation  with  the  corresponding  values  obtained  by  using  Box's  asymptotic 

-13 

expression  of  order  n  .  Upper  percentage  points  of  the  distribution  of  >3 
for  a  =  0.01,  0.05  and  s  *  3(1)6  are  also  given  by  Nagarsenker  and  Das  (1975) 
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Table  3 


Comparison  of  the  Pearson  Type  Approximation  with  the  Asymptotic  Expression 
for  the  Distribution  of 


s  =  5 _ 

s  =  S 

n 

*^3 

“l 

“2 

^^3 

“l 

“2 

15 

2.763 

0.05 

0.0496 

- 

- 

- 

15 

3.265 

0.01 

0.0099 

- 

- 

- 

21 

1.895 

0.05 

0.0498 

4.565 

0.05 

0.0490 

21 

2.238 

0.01 

0.0100 

5.093 

0.01 

0.0097 

41 

0,928 

0.05 

0.0500 

2.160 

0.05 

0.0499 

41 

1.095 

0.01 

0.0100 

2.409 

0.01 

0.0100 

51 

0.739 

0.05 

0.0500 

1.711 

0.05 

0.0499 

51 

0.873 

0.01 

0.0100 

1,908 

0.01 

0.0100 

In  the  above  table,  is  the  value  of  a  obtained  by  using  the  Pearson  type 
approximation  whereas  ^2  is  ^he  value  of  a  obtained  by  using  Box’s  asymptotic 
expansion  of  order  n'^^  This  table  indicates  that  the  accuracy  of  Pearson 
type  approximation  is  sufficient  for  practical  purposes. 
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6.  Test  Specifying  the  Covariance  Matrix 

The  modified  likelihood  ratio  statistic  for  testing  the  hypothesis 
E  =  is  given  by 

-  (e/n)^"  etrC-AI^^).  (6.1) 

The  modified  likelihood  ratio  test  statistic  is  obtained  from  the  likeli¬ 
hood  ratio  test  statistic  by  changing  N  to  n.  The  moments  of  are  seen 
to  be 

E(X^)  =  (e/n)®^"  |l  1"^  ll  +  hi 
4  'o'  o' 

s 

X  n  {r(n+nh+l-l)/r(n+l-i)}.  (6.2) 

i=l 

Anderson  (1958)  derived  the  likelihood  ratio  statistic  for  testing  the 
hypothesis  that  the  covariance  matrix  is  equal  to  a  specified  matrix  when 
the  underlying  distribution  is  real  multivariate  normal.  The  distribution 
of  can  be  approximated  with  a  Pearson  Type  I  distribution  using 

the  first  four  moments,  where  b  is  a  suitably  chosen  integer.  Using  the 
above  approximation,  Krlshnalah,  Lee  and  Chang  (1975,  1976)  computed  the 
percentage  points  of  the  distribution  of  X^  =  -2  log  for  s  =  2(1)10, 
a  =  0.01,  0.05,  M=l(l)20(2)30,  where  M  -  n-s-1  and  P[X^<c^1h^]  -  (1-a). 

These  percentage  points  are  reproduced  in  Table  10. 
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7.  Test  for  Multiple  Homogeneity  of  the  Covariance  Matrices 

Lst  ^l’****^q  independently  distributed  as  complex  p-variate  normal 

with  mean  vectors  and  covariance  matrices  E,, . E  ,  respectively. 

--*■  -4  11  qcj  ^  ^ 

Also,  let  (j  •  1 , . . .  ,N^)  be  the  independent  observation  on  In  this  sec¬ 

tion,  we  study  the  Pearson  type  approximation  to  the  distribution  of  the 


likelihood  ratio  statistic  for  testing  where 


E,  »...»E  „ 

11 


^q^+l.q^+l“"-’^q*.q* 


qq 


qQ  *  0*  q^^  “  qj«  q^  *  q  q^  “  q^'*’»**'*'qj*  The  modified  likelihood  ratio 
statistic  (obtained  by  changing  to  n^  in  the  likelihood  ratio  test  statis¬ 
tic)  for  testing  is  given  by 


q  ,  n. 

n  |A  /n  I 
i*l 


A  /n.  I 

j=l  i*q*  ,+l  ^ 

J-1 


!  i. 


(7.1) 


where 


*  ‘’J 

"i  =  ‘  I  "i 


and 


l-q*_^+l 


^ii  ’  ^  ^!ij  -  ?i.^^?ij  '  •  ?i. 


j*l 

The  moments  of  are  given  by 


j-1' 
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1, 


ECXg) 


*  “I 


q  I 
n  (n^) 

a*l 


P  d 

n  n 

i-l  a-1 


r(n  +hn  +l-i) 

_ S _ — S _ 

r(n  +l-i) 

6 


r(n  +l-i) 

_ g 

r(n  +hn  +l-i) 

a  g 


Using  the  first  four  moments  of  X^,  the  distribution  of  X^  can  be  approxi¬ 
mated  with  a  Pearson's  Type  I  distribution  where  b  Is  a  suitably  chosen 
Integer.  This  approximation  was  used  by  Krlshnalah,  Lee  and  Chang  (1975, 
1976)  to  compute  approximate  percentage  points  of  the  distribution  of 
X^  -  -2  log  X^  for  n^  =  n^,  q  »  dk  (l.e.,  there  are  k  populations  In 
each  of  the  d  groups).  These  points  are  reproduced  In  Table  11  for  d=l. 

In  Tables  4  and  5,  we  compare  the  values  obtained  by  the  Pearson  Type 
approximation  for  the  distribution  function  of  with  the  corresponding 

values  obtained  by  using  the  Box's  asymptotic  expansion  up  to  terms  of  order 
-13 

n  .  In  these  tables,  the  constant  c^  Is  defined  as 

P[-2  log  Xj^Cj  =  (l-g).  (7.3) 


Also,  Is  the  value  of  g  If  we  use  the  Pearson  type  approximation  whereas 

-13 

is  the  value  of  a  if  we  use  the  asymptotic  expression  of  order  n 
Tables  4  and  5  indicate  that  the  accuracy  of  the  Pearson  type  approximation 
is  sufficient  for  practical  purposes. 

In  the  real  case,  Wilks  (1932)  derived  the  likelihood  ratio  statistic 
for  testing  the  homogeneity  of  the  covariance  matrices  whereas  Krishnaiah 
and  Lee  (1976)  discussed  how  certain  tests  of  hypotheses  on  linear  struc¬ 
ture  of  the  covariance  matrices  can  be  reduced  to  the  problem  of  testing  for 
the  multiple  homogeneity  of  the  covariance  matrices. 
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Table  4 


Comparison  of  the  Pearson  Type  Approximation  with  the  Asymptotic  Expansion 
for  the  Distribution  of  when  d  »  1 


_ P  =  3 _ 

_ P  °  _ _ 

"o 

q 

^=5 

“1 

“2 

*=5 

“1 

“2 

10 

2 

19.82 

0.05 

0.0501 

33.00 

0.05 

0.0502 

10 

2 

25.40 

0.0100 

40.21 

0.01 

0.0101 

10 

6 

69.68 

0.05 

0.0500 

12.14 

0.05 

0.0492 

10 

6 

79.11 

0.0100 

13.39 

0.01 

0.0098 

15 

2 

18.72 

0.05 

0.0501 

30.33 

0.05 

0.0500 

15 

2 

23.99 

0.01 

0.0100 

36.93 

0.01 

0.0100 

15 

6 

66.70 

0.05 

0.0500 

113.77 

0.05 

0.0498 

15 

6 

75.69 

0.01 

0.0100 

125.48 

0.01 

0.0099 

20 

2 

18.23 

0.05 

0.0500 

29.18 

0,05 

0.0500 

{20 

2 

23.35 

0.01 

0.0100 

35.52 

0.01 

0.0100 

•20 

j 

6 

65.34 

0.05 

0.0499 

110.45 

0.05 

0.0499 

I20 

6 

74.13 

0.01 

0.0100 

121.78 

0.01 

0.0100 
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8.  Simultaneous  Tests  for  the  Homogeneity  of  Populations 
In  this  section,  we  discuss  the  likelihood  ratio  test  for  the  homo¬ 
geneity  of  complex  multivariate  normal  populations.  The  hypothesis  of  the 
homogeneity  of  the  q  complex  multivariate  distributions  defined  in  Section 
7  is  equivalent  to  the  hypothesis  where 


(8.1) 


The  modified  likelihood  ratio  statistic  for  testing  is  given  by 


nP"  n  IG  1  ^ 

i=l 


6  q  Pn,  q  _ 

n  n,  ]G+l  N  (z  -z..)  (z  -z..)'  I  " 

1-L  ^  1-1  ^  ' 


where  n  *  f  n  ,  n  =  N.-l, 
i=l  1 


M  1. 

N  -  N  .z..  «  i  I  I  z  ,  Zj  [  z 

i=l  1  '  1  j  -iJ  j=i  -iJ 


(8.2) 


1  =  I  (!ir?i.^(!ir!i.^‘ 


j=i 


G  -  y  G. 

L  ,  . 

i=l 


The  moments  of  are  given  by 


.  q  phn.  p  q  r(n,+hn.+l-i) 


r (n+q-i) 


(8.3) 


The  distribution  of  is  approximated  with  a  Pearson  Type  1  distribution 

where  b  is  a  suitably  chosen  integer.  Using  this  approximation,  percentage 
points  of  the  distribution  of  X^  =  -2  log  X^  are  computed  by  Chang,  Krishnaiah 
and  kee  (1975)  for  a  »  0.01,  0.025,  0.05,  O.lO;  n^  •  **0 

p»  1,  2,  3,  4  and  M  »  n^-P»  1(1)  20,  25,  30.  Table  12  gives  the  value  of  Cg 
for  a  •  0.05,  0.01  where  Cg  is  given  by 

17 


4 


(1-a). 


(8.4) 


To  check,  for  the  accuracy  of  the  entries  in  Table  12,  we  compared 

some  of  these  values  with  the  corresponding  values  obtained  by  using  Box's 

-13 

asymptotic  series  of  order  n  .  These  comparisons  are  given  in  Table  6. 


Table  6 

Comparison  of  the  Pearson  Type  Approximation  with  the  Asymptotic  Expansion 


_ 1  _ 

L _ ' _ 

n 

0 

q 

r — 

<=6 

“1 

“2 

1 

^=6 

“1 

“2 

7 

3  i 

28.03 

0.05 

i 

0.0499 

50.49 

0.05 

0.0497 

7 

3 

34.15 

O.Ql 

0.0100 

58.80 

0.01 

0.0098 

10 

3 

27.45 

0.05 

0.0500 

48.06 

0.05 

0.0498 

10 

3 

33.42 

0.01 

0.0100 

55.92 

0.01 

0.0099 

15 

3 

27.04 

0.05 

0.0500 

46.46 

0.05 

0.0499 

15 

3 

32,90 

0.01 

O.OlOO 

54.04 

0.01 

0.0100 

20 

3 

I  26.84 

0.05 

0.0500 

45.73 

0.05 

0.0500 

20 

3 

32.67 

0.01 

0.0100 

I  53.18 

1 

0.01 

0.0100 

1/b 

In  the  above  table,  is  the  value  of  a  obtained  by  approximating  with 

Pearson's  type  I  distribution  whereas  02  is  the  value  of  a  detained  by  using 
Box's  asymptotic  series.  Table  6  indicates  that  the  accuracy  of  the  values 
of  is  good. 
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9.  Test  Specifying  the  Values  of  the  Covariance  Matrix  and  Mean  Vector 
In  this  section,  we  consider  the  distribution  of  the  likelihood 
ratio  statistic  for  testing  the  hypothesis  where 


'h  “  i!o 


and  and  are  known.  The  likelihood  ratio  statistic  for  testing 
is  given  by 


pN,  ,  N 

X7  =  (e/Nj)  |G^  Z-^\ 


etr[-r^  {G^+N^(z^^-U^)(zi  -Uo)'}]  (9.1) 


The  moments  of  are  given  by 

.  phN,  ,  P  r(N  -i+Nh) 

=  <ir)  ^ 

^  (l+h)  ^ 

Using  the  first  four  moments,  Chang,  Krishnaiah  and  Lee  (1975,  1977) 
approximated  the  distribution  of  with  Pearson  Type  I  distribution 

where  b  is  a  suitably  chosen  integer.  This  approximation  is  used  to 
compute  the  values  of  c^  for  M=*nj-p-l  «  1(1)20(2)30,  njS'Nj-l  and  p=2,3,4,5,6 
where 


P[X^  <  c^lH^]  =  (1-a). 


and  \y  =  -2  log  X^.  These  values  are  given  in  Table  13  for  a  =  0.05,  0.01, 
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10,  Applications  in  Time  Series  in  the  Frequency  Domain 

In  this  section,  we  discuss  as  to  how  the  Likelihood  ratio  test 
procedures  on  the  covariance  matrices  of  the  complex  multivariate  normal 
populations  can  be  used  in  the  area  of  Inference  on  multiple  time  series. 

Let  X'(t)  *  (X^(t),  X^(t))  (t  ■  1,...,  T)  form  a  Gaussian, 

stationary,  multiple  time  series  with  zero  means  and  covariance  matrix 
R(s)  =  (Rjj^(s))  where  R^j^(s)  =  E{Xj(t)  X^(t  +  s)}  and  Xj(t)  is  of  order 
pj  1.  The  spectral  density  matrix  of  the  above  time  series  is  given  by 
F(iii)  =  (F^j(ui))  where 


-f,'  I  *  *"' 


(10.1) 


A  well  known  estimate  (e.g.  see  Brillinger  (1974)  of  F(u)  is  F(w)  = 


(Fi,.(w))  where 
*■1 


(10.2) 


and  ffl  is  a  suitably  chosen  integer,  ta  Eq.  (10.2) 

1  T 

°  - -  I  X  (t)  exp  (-itX), 

/27tT  t=l  ' 

and 

lij(A)  =  Z^(l)  Z’(A). 

Goodman  (1963b)  and  Wahba  (1968)  showed  that  A(uj)  =  (2m+l)  F(uj)  is  approxi¬ 
mately  distributed  as  complex  Wishart  distribution  with  (Im+l)  degrees  of 
freedom  and  E(F(ai))  *  Ffw). 

Now,  let  H2(w),  H2(«)  and  H^(a))  denote  the  following  hypotheses: 
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(10.3) 


H2{w):  Fjjj(u))  =  0  (i  j  =  1,  ....  q) 
F((ij)  =  o^F^Cu) 

H^(oj):  F((jj)  =*  F^Coj) 


where  a  Is  uoknotm  and  F^(a))  is  known.  Let  the  statistics  12(1^)),  L^Cto), 


and'L^(oj)  be  defined  as  follows; 


L2(,)  =  JAMJ - 

n  |a^^(w)| 


1=1 


(10.4) 


L-(tj)  = 


|A(oj)  F"^(aj)| 


{tr  A(u))  F  ^(w)/s}^ 
0 


(10.5) 


L^(u)  =  (e/(2nri-l)j®^^®*'^^  1a(uj)  F'^(.^)| 


-1,  ,|(2m+l) 


(10.6) 


X  etr  (-A(u))  F^^(u!)). 


where  (2nH-l)A^j  (oj)  =  Also,  let  L^(u,)  =  -2  log  L^(u)  for  1  -  2,3,4. 


■ij 


The  hypothesis  H2(w)  is  accepted  or  rejected  according  as 


12(01)  >  d; 


where 


P[L2(iij)  <  d2|H2(ui)]  =  (1  -  a). 

We  accept  or  reject  H2(o))  according  as 


12(01)  >  dj 


(10.7) 


(10.8) 


where 


P[L2(oi)  <  d2|H2(u.)]  -  (1  -  a). 


(10.9) 
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Similarly,  the  hypothesis  is  accepted  or  rejected  according  as 

L, (a)  >  d, 

H  4 

where 

PtL^(a)  <  d^lH^Ca)]  =  (1  -  a).  (10.10) 

Since  A(a)  is  approximately  distributed  as  the  complex  Wishart  matrix, 

approximate  values  of  d^,  d^,  and  d^  can  be  obtained  from  Table  8,  Table  9, 

and  Table  10,  respectively. 

k  k  k 

Next,  let  H  =  '^  H  (a.),  H  -O  H-(a  )  and  H  -  O  H  (a  ) 
j  =  l  ^  J  j  =  l  ^  J  ^  1 

where  ap...,aj^  are  widely  separated.  Then,  we  accept  or  reject  H2 

according  as 

T2  >  (10.11) 

where 

P[T2  <  d3lH2]  =  (1-a)  (10.12) 

and  T2  *  max  (L2 (uj) , . . . ,L2 («j^) )  .  An  alternative  procedure  is  to  accept  or 
reject  H2  according  as 

T*  $  d^  (10.13) 

where 

P[T2  <  °  (10.14) 

and 

k  . 

T*  =  n  L  (a  ). 
j-1  ^ 

Since  a^,...,aj^  are  widely  separated,  L2(aj)...,  L2(a^)  are  distributed 

Independently.  So  the  critical  values  d^  and  d^  can  be  computed  by  using  the 

methods  discussed  in  this  chapter.  We  can  propose  similar  procedures  to 

test  H_  and  H,  . 

3  4 
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Let  Hj  F(a)j^) 


=  ...  =  F(u)j^),  where  the  frequencies  are 

"sufficiently"  wide  apart.  Also,  let  L^Cui)  =  -2  log  ,  where 

I  lF<a>  )l(2nri-l) 

L  _ 

'  II  F(.  l/kl'^^^nr.l) 

1=1 

Then,  we  accept  or  reject  accordingly  as  f  d*  where 
P[L5  <  dllHj]  =  (1-a). 

Since  , . . .  ,ojj^  are  "sufficiently"  wide  apart,  F(a)j^),...,  F(a)j^)  are  dis¬ 
tributed  independently.  Also,  (2id+1)  F(w^)  is  distributed  approximately 
as  the  complex  Wishart  matrix  with  (2m+l)  degrees  of  freedom  for  i=l,...,k. 
Hence,  the  values  of  d*  can  be  obtained  from  Table  11. 

Next,  let  X^(t)  =  (t) , . . . ,X’^ (t) ) ,  (t=l . T^)  be 

a  Gaussian,  stationary,  multiple  time  series  with  zero  means  and  covariance 
matrix  Rj^(s)  and  spectral  density  matrix  F^(u),  where  R^Cs)  =  ^'^iuv^®^^ 

^uv(")  =  x;^(t+s)}. 

Also,  let  Xj^(t)  , . . .  ,Xj^(t)  be  distributed  independently  and  Xj^(t) 
be  of  order  pxl  for  i=l,...,k.  Let  the  estimate  F^(a.)  of  F(^)  be  defined 
in  a  similar  way  as  F(w).  Here,  (2m^+l)  Fj^(a))  is  distributed  approximately 
as  the  complex  VJishart  matrix  with  2m^+l  degrees  freedom.  The  hypothesis 
H,(w);  F.  (uj)  =  ...  =  F,  (u>)  is  tested  as  follows.  We  accept  or  reject  H. 

o  t  K  g 

accordingly  as 

Lg(i.))  >  dg 
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where 


P[L^(^)  <  d*|Hgl  =  (1-a), 
=  “2  log  L^(a)), 


L^(u) 


k  . 

n  |f  (uj) 

1=1 


(2m^+l) 


k  .  m 

I  (2m^+l)F.(a,)/m^|  ° 


k 

and  m  =  2(  I  m. )  +  k.  The  critical  values  dc  can  be  obtained  from  Table 
°  1=1  ^  ^ 

11  when  the  m^'s  are  equal. 

We  will  now  illustrate  the  usefulness  of  some  of  the  tables  in  this 
paper  with  vibration  data  on  a  C-5A  transport  aircraft. 

Vibration  measurements  have  been  taken  on  the  cargo  deck  of  a  C-5A 
transport  aircraft  to  provide  information  about  the  dynamic  environments 
that  cargo  must  survive  in  transit  and  to  understand  better  the  distribu¬ 
tion  and  transmission  of  vibrational  energy  throughout  the  aircraft  struc¬ 
ture.  Measurements  have  been  taken  over  certain  periods  by  locating 


accelerometers 

at  different  locations 

on  Che  cargo  deck. 

We 

will  treat 

each  location 

as  a  variable.  Data  on 

the  following  variables 

were  taken: 

Longitudinal • 

Lateral 

Directional 

Variables 

Location 

Location 

Orientation 

1  (FRV) 

Forward 

Right 

Vertical 

2  (FRL) 

Forward 

Right 

Lateral 

3  (FLV) 

Forward 

Left 

Vertical 

4  (FLL) 

Forward 

Left 

Lateral 

5  (ARV) 

Aft 

Right 

Vertical 

6  (ARL) 

Aft 

Right 

Lateral 

7  (ALV) 

Aft 

Left 

Vertical 

8  (ALL) 

Aft 

Left 

Lateral 
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The  basic  unit  of  measurement  is  the  acceleration  due  to  gravity  (Ig. = 

2 

980  cm/sec  ).  Let  the  spectral <fensity  of  the  data  on  the  above  8  variables 
at  frequency  u  be  denoted  by  F(a))  and  let  the  corresponding  population 
spectral  density  matrix  be  denoted  by  F(<o).  The  sample  spectral  density 
matrix  at  frequency  is  F(u)j)  -  +  iS^^  where  =  0.15907HZ, 

a  “  2, 3, 4, 5  and 
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Let  F(a>)  be  partitioned  as 


F(u) 


ll(<.) 

^2(“> 

Fi3(.) 

F,4(.) 

21(0.) 

F22(<^) 

F23(.) 

F2,(a.) 

31  (^) 

F32(<.) 

F33(.) 

F34(.) 

4i<“) 

F^2(“) 

"43<“> 

of  order  2x2 

for  i  *  1 

,2,3,4. 

=  -2  log 


L2(Wj)  where 


(2nri-l)F(a)j)| 


n  |(2nri-l)F^^(u  )| 
1=1  ^ 


and  (2in+l)  =  19.  The  values  of  L  (w^  in  this  case  are  found  to  be  47.760, 
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31,667,  33.684,  35.646  and  37.738  respectively.  The  value  of  the  critical 
value  d2  for  n=19,  s=8  and  q=4  from  Table  8  is  found  to  be  4.217  at  5% 
significance  level.  Since  the  computed  values  of  are  greater  than 

the  value  from  the  table,  we  conclude  that  the  sets  (1,2),  (3,4),  (5,6), 
(7,8)  of  variables  are  not  independent  for  each  of  the  five  frequencies 
considered. 

Next,  we  computed  the  value  of  12(01^)  =  -2  log  L2(u)j) 

where 


|(2nri-l)F(a)  )1 

L2(ui  )  =  .  ^ 

{  (2m+l)tr  F(a)j)/s} 

s»9,  (2nH-l)  =  19.  It  is  found  that  the  values  of  Lj(u)j),  in  this  case, 
are  78.321,  65.267,  63,226.  59.019  and  61.852  respectively.  The  critical 
value  d^  for  n»19  and  s»8  is  found  to  be  5.142  at  5%  level.  So,  we 
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reject  (individually)  the  hypotheses  that  “  o  Ip  j  “  1»2,3,4,5. 

For  the  hypothesis  we  will  consider  only  the  first  four  variables. 

Let  spectral  density  matrix  of  the  data  on  the  first  four  variables  at 
frequency  be  devoted  by  F(u)^)  and  let  the  corresponding  population  spectral 
density  matrix  be  denoted  by  F(aj^).  The  sample  spectral  density  matrices  at 
frequencies  0.15907  Hz,  0.47721  Hz  and  0.79535  Hz  are  F(0. 15907)  =  +  i  S, 
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i  and 

F(0. 79535) 

“  "30  ^  ^ 

res 

"io  = 

i2O4*at0 

.0027030 
•  001«>000 
•  002(»990 

••0000933 

>s;o>ioos 
.•02(440 
. 0316000 
.0036660 

.6033170 

-.0000933 

.0036666 

•6606696 

®11  “ 

6,<t;0aOQtt 
• JOMkOOO 

•  «aoO«iQ39 

••0044000 

0*0000000 

•«00497$0 

••0003703 

'.0106100 
. 0064140 

0. 0040000 
'.0001100 

•6604039 

•0003763 

•0007700 

9.0600000 

^20  “ 

•oaoitis 

«9usr«a 

•-«<ifl09tor 

■•V 

•0000011 

•0007030 
•  0009133 
••0000139 

.0104100 

.0004133 

.0346000 

•.0006346 

•*1005107 

••6006139 

•*6006396 

•000S513 

"21  • 

0*0000000 
•OOObO^S 
•  •00t22<)0 
•«00««|000 

••0000049 

••ooooooo 

••0002020 

••0000029 

.0012140 
. 0002020 

0. OOOOOOO 
'.0001446 

•6004696 

•0000029 

.0001996 

0.0000000 

"30  • 

••]/7990 

••0002297 

«0I9%000 

•0007)%« 

••0002297 

•0000092 

••0003207 

••0000763 

. 0146000 
'.0003201 
.0214400 
. 0010120 

*0007346 

•*0006763 

*0030720 

•0007642 

"31  = 

OeOOOOOOO 

«900<i««0 

••0002979 

•••aa%923 

••0004046 

9*0000000 

••0003417 

•0000041 

.0002414 

.0003611 

0. OOOOOOO 
'.0003601 

*0004523 

-*0000043 

•0003403 

0*0000009 

29 


J 

] 

] 

] 


J 


10^  =  0.15907  Hz,  -  0.47721  Hz,  =  0.79535  Hz,  (2nri-l)  =  19.  The  value  of 
L c  is  found  to  be  224.72.  The  critical  value  d*  for  n  =19,  p=4  and  k=3  from 
Table  li  is  found  to  be  50.93  at  the  5%  significance  level.  Since  L ^  is 
greater  than  the  value  from  the  table,  we  conclude  that  the  spectral 
density  matrices  F(ojj^),  F(w2)  and  FCui^)  are  significantly  different  from 
each  other. 
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1. 


Appendix 


Tables  7  through  13  give  percentage  points  of  various  distributions 
discussed  in  this  chapter.  These  tables  are  useful  in  implementation 
of  the  likelihood  ratio  test  procedures.  A  description  of  these  tables 
is  given  below. 
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Table  7:  Percentage  Points  of  the  Distribution  of  the  Determinant  of 
the  Complex  Multivariate  Beta  Matrix 

The  entries  in  this  table  give  the  values  of  where 
P[Cj  <  Cj]  -  (1  -  a), 

=  -{ (2n+q-p)log  and  U  =  | (Aj+A2)"^ I . 

Here  Aj^  and  A2  are  distributed  independently  as  central  pxp  complex 
Wishart  matrices  with  n  and  q  degrees  of  freedom  respectively.  These 
entries  in  this  table  are  reproduced  from  a  technical  report  by  Lee, 
Krishnaiah  and  Chang  (1975). 


Here,  A  =  :  pq*pq  is  the  central  complex  Wishart  matrix  with  n 

degrees  of  freedom  and  E(A)  =  nT,  and  T  “  Also,  A^^  and  are  of  order 

pxp.  In  addition,  =  0  for  i^* j=l , 2 , . . .  ,q  ,  M=n-s-3  and  s=pq. 
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Krlshnalah,  Lee  and  Chang  (1976)  gave  the  values  of  for 
a  “  0,05,  M  =  1(1)4(2)16,20,24,30,  p=l,2,3,  and  q-3,4,5.  These  entries 
are  Included  in  Table  8  with  the  kind  permission  of  Biometrika  Trustees. 
The  remaining  entries  in  the  table  are  reproduced  from  a  technical  report 
by  Krishnaiah,  Lee  and  Chang  (1975). 

Table  9;  Percentage  Points  of  the  Likelihood  Ratio  Test  Statistic 
for  Sphericity 

The  entries  in  this  table  are  the  values  of  where 

P(X3  1  C3]  =  (1  -  Qt), 

Xj  -  -2  log{  lAl“^l/(trAI~Vs)®}, 

and  M=n-s-3.  Also,  A:sxs  is  distributed  as  the  central  complex  Wishart 

2  2 

matrix  with  n  degrees  of  freedom  and  E(A)  «  no  E  where  0  is  unknown  and 

o 

is  known.  Upper  5%  points  of  the  distribution  of  X^  are  given  in 
Krishnaiah,  Lee  and  Chang  (1976)  for  s=7(l)10  and  M=1 ( 1) 5 , 7 , 10 , 15 , 20 , 
30(5)50,60.  These  entries  are  reproduced  in  Table  9  with  the  kind 
permission  of  Biometrika  Trustees.  The  remaining  entries  in  the  table 
are  reproduced  from  a  technical  report  by  Krishnaiah,  Lee  and  Chang  (1975). 

Table  10:  Percentage  points  of  the  Likelihood  Ratio  Statistic  for 
Specifying  the  Covariance  Matrix 

The  entries  in  this  table  give  the  values  of  c  where 

4 

PtX^  £  c^]  -  (1  -  a), 

X.  »  (e/n)^"  !aZ  ^1”  etr(-AZ  ^)  and  \  «  -2  log  X,. 

4  O  O  4  4 
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In  the  above  equation,  A:sxs  is  distributed  as  the  central  Wishart  matrix 

with  n  degrees  of  freedom  and  E(A)  =  nZ  where  Z  is  known.  Krishnaiah, 

o  o 

Lee  and  Chang  (1976)  gave  upper  5%  points  of  the  distribution  of  for 
s=2(l)7,  M=1 (1) 5, 7 , 10 , 15 ,20,30  where  M“n-s-l.  These  percentage  points 
are  reproduced  in  Table  10  with  the  kind  permission  of  the  Biometrika 
Turstees.  The  remaining  entries  in  the  table  are  reproduced  from  a 
technical  report  by  Krishnaiah,  Lee  and  Chang  (1975), 


Table  11.  Percentage  Points  of  the  Likelihood  Ratio  Tests  Statistic 
for  the  Homogeneity  of  the  Covariance  Matrices 


The  entries  in  this  table  are  the  values  of  c^  where 
Ptx's  <  c^]  =  (1  -  a). 


where  ^ii»****Aqq  3re  distributed  independently  as  central  pxp  complex 
Wishart  matrices  with  nj,..,,n^  degrees  of  freedom  respectively, 

E(Aj ^/nj)=. . .=E(A^j^/n^)  and  n=n|+...+nq.  Upper  5%  points  of  the  distri¬ 
bution  of  X^  are  given  in  Krishnaiah,  Lee  and  Chang  (1976)  for  p  =  3,4;  q 
2(1)6, 8;  n^=5(l)20,25,30  where  nj=.  ♦ -“nq^np*  These  points  are  reproduced 
in  Table  11  with  the  kind  permission  of  Biometrika  Trustees.  The 
remaining  entries  in  this  table  are  reproduced  from  a  technical  report 
by  Krishnaiah,  Lee  and  Chang  (1975). 
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Table  12:  Percentage  Points  of  the  Likelihood  Ratio  Test  Statistic 

for  the  Homogeneity  of  Complex  Multivariate  Normal  Populations 

The  entries  in  this  table  are  the  values  of  c  where 

o 

P[Xg  _<  CglHg]  -  (1  -  a), 

Xg  *  -2  log  Xg  and  X^  is  the  likelihood  ratio  statistic  for  testing  the 
hypothesis  of  the  homogeneity  of  q  complex  p-varlate  normal  populations. 
The  statistic  X^  Is  defined  In  Section  8.  In  the  table,  q  denotes  the 
number  of  populations  M~N^-p  and  Is  the  cosason  size  of  various  groups. 
Chang,  Krlshnalah  and  Lee  (1977)  gave  upper  5%  points  of  the  distribution 
of  X^  for  q«2(l)5  and  M”l(l)20,25,30.  These  percentage  points  are  repro¬ 
duced  In  Table  12  with  the  kind  permission  of  the  North-Holland  Publishing 
Company.  The  remaining  entries  in  this  table  are  reproduced  from  the 
technical  report  by  Chang,  Krlshnalah  and  Lee  (1975). 


Table  13;  Percentage  Points  of  the  Likelihood  Ratio  Test  for  Z  *  Z 

and  u  >■  u  . 

.  -O 


The  entires  in  this  table  are  the  values  of  e ^  where 

PIX;  <  Cj^lH^l  -  (1  -  a),  X^  -  -2  log  X^ 
and  X^  Is  the  likelihood  ratio  statistic  for  testing  the  hypothesis 
Hy  where 


”7=^1  ■  ^o’  “1  "  “o* 

Here  and  are  respectively  the  mean  vector  and  covariance  matrix 

of  p-varlate  complex  normal  distribution  and  u  and  Z  are  known.  Also, 

.00  ’ 

M»N^-p-2  and  N^  is  the  sample  size.  Chang,  Krlshnalah  and  Lee  (1977) 
gave  the  values  of  Cj  for  a  •  0.05,  M  »  1(1)20(2)30  and  p  »  2(1)6. 

These  values  are  reproduced  in  Table  13  with  the  kind  permission  of  North- 
Holland  Publishing  Company.  The  remaining  entries  are  reproduced  from 
the  technical  report  by  Chang,  Krlshnalah  and  Lee  (1975). 
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Some  of  the  entries  in  Tables  2-5  are  reproduced  from  Krishnaiah , 

Lee  and  Chang  (1976)  with  the  kind  permission  of  the  Biometrlka  Trustees. 
Whereas  some  of  the  entries  in  Table  6  are  reproduced  from  Chang,  Krish¬ 
naiah  and  Lee  (1977)  with  the  kind  permission  of  North-Holland  Publishing 
Company. 
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Percentage  Points  of  the  Distribution  of  the  Determinant  of  the  Complex  Multivariate  Beta 
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Table  7  (continued) 
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Table  7  (continued) 
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Table  7  (continued) 
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